DARK HALOS: 
THE FLATTENING OF THE DENSITY CUSP BY 
DYNAMICAL FRICTION 
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ABSTRACT 

N-body simulations and analytical calculations of the gravitational collapse in an expanding 
universe predict that halos should form with a diverging inner density profile, the cusp. There are 
some observational indications that the dark matter distribution in galaxies might be character- 
ized by a finite core. This 'core catastrophe' has prompted a search for alternatives to the CDM 
cosmogony. It is shown here that the discrepancy between theory and observations can be very 
naturally resolved within the standard CDM model, provided that gas is not initially smoothly 
distributed in the dark matter halo, but rather is concentrated in clumps of mass > 0.01% the 
total mass of the system. Dynamical friction acting on these lumps moving in the background 
of the dark matter particles, dissipates the clumps orbital energy and deposits it in the dark 
matter. Using Monte-Carlo simulations, it is shown that the dynamical friction provides a strong 
enough drag, and that with realistic baryonic mass fractions, the available orbital energy of the 
clumps is sufficient to heat the halo and turn the primordial cusp into a finite, non-diverging core 
— overcoming the competing effect of adiabatic contraction due to gravitational influence of the 
shrinking baryonic component. Depending on the initial conditions, the total density distribution 
may either become more or less centrally concentrated. Possible consequences of the proposed 
mechanism for other problems in the CDM model and for the formation and early evolution of 
the baryonic component of galaxies are also briefly discussed. 



Subject headings: galaxies: evolution - galaxies: ISM - galaxies: kinematics 
structure - hydrodynamics 



dynamics - galaxies: 



1. Introduction and motivation 

The formation and structure of dark matter 
(DM) halos around galaxies and clusters of galax- 
ies poses one of the great challenges to theories 
of structure formation in an expanding universe. 
The basic theoretical paradigm was set by Gunn & 
Gott (1972) and Gunn (1977), and a more detailed 
analysis considering the dependence on the pri- 



mordial perturbation field was done by Hoffman & 
Shaham (1985). The first N-body simulations de- 
signed to probe the structure of dark halos indeed 
confirmed the basic findings of Hoffman and Sha- 
ham (Quinn, Salmon & Zurek 1986; Frenk et al. 
1988), but could not resolve the innermost struc- 
ture of the dark halos. Higher-resolution simula- 
tions showed evidence of a cuspy density profile of 
the form p(r) oc r^" with a > 1, down to the res- 
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olution limit of the simulations (Dubinski & Carl- 
berg 1991; Warren et al. 1992; Crone et al. 1994), 
and further limited the exponent of the innermost 
density profile to a k, (1 — 1.5) (Klypin et al. 2000; 
Navarro, Frenk & White 1997, hereafter NFW; 
Cole & Lacey 1996; Tormen, Bouchet & White 
1997; Huss, Jain & Steinmetz 1999; Fukushige & 
Makino 1997; Moore et al. 1999; Jing & Suto 
2000). NFW suggested a particular parametric 
density profile, claimed to be a universal one (in- 
dependent of halo mass), and characterized by an 
a = 1 cusp and r~^-like decay at large radii. On 
the theoretical side it has been demonstrated that 
a cuspy density profile arises inherently from the 
cold gravitational collapse in an expanding uni- 
verse (Lokas & Hoffman 2000). 

Whether the predicted density cusp is sup- 
ported or ruled out by observations is subject 
to debate. Moore (1994) and Flores & Primack 
(1994) argued that the observed HI rotation curves 
of gas-rich dwarf galaxies favor a core-halo struc- 
ture rather than a diverging inner density pro- 
file. This has been confirmed by more recent ob- 
servations of De Blok et al. (2001). The DM 
distribution inferred for spiral galaxies also show 
large constant density cores (e.g., Salucci & Burk- 
ert 2000). Moore et al. (1999) coined this dis- 
agreement as the 'core catastrophe,' using it as 
a major argument against the cold dark matter 
(CDM) model. Indeed, alternative theories for 
structure formation have been proposed, e.g. self- 
interacting dark matter (SIDM: see, e.g.,Spergel & 
Steinhardt 2000). Nevertheless, it has been sug- 
gested that the inconsistency may reflect the finite 
resolution of the observations that has not been 
properly accounted for in the analysis of the HI 
rotation curves (van den Bosch et al. 2000: van 
den Bosch & Swaters 2000). Moving on to clus- 
ters of galaxies, the cuspy density profile agrees 
with the recent determination of the mass profile 
from ASCA X-ray observation of A2199 and A496 
(Markevitch et al. 1999). These cannot be ac- 
counted for in the SIDM scenario (Yoshida et al. 
2000). In addition, gravitational systems are ulti- 
mately unstable to heat transport and SIDM mod- 
els will eventually undergo gravothermal catastro- 
phe, resulting in even more centrally concentrated 
structures: (Moore et al. 2000; Kochanek & White 
2000). 

Assuming that the theory of cold collapse in an 



expanding universe is correct, that halos form with 
a density cusp, and that observations of galaxies 
indicate that the present day structure has a finite 
non-diverging core, we discuss a mechanism for 
resolving the above conflict. Close inspection of 
the problem at hand reveals that adding a mech- 
anism that enables 'heat' transport to the inner 
part of the halo would resolve the problem — de- 
stroying the inner cold region where the "temper- 
ature" inversion occurs. This will lead to the puff- 
ing up of the central regions and to the flattening 
of the density cusp. The aim of this paper is to 
show that such a mechanism can be naturally for- 
mulated within the CDM cosmogony without ap- 
pealing to 'new physics. 'In addition, as opposed to 
SIDM we have a two component system - with 
one component (the halo) that is always expand- 
ing and the other (baryonic) component that is 
always shrinking. The baryonic component be- 
coming centrally concentrated is not an embar- 
rassment, since galaxies do have disk-bulge-central 
mass components that arc more centrally concen- 
trated than the observed halo structures. In ad- 
dition, part of the baryonic mass can be lost in 
subsequent outflow (which is efficient in precisely 
those galaxies without highly concentrated bary- 
onic components). 

The basic scenario envisaged here is one where 
the gas {i.e. baryon) distribution becomes lumpy 
because of gravitational collapse and cooling on 
its Jeans scale. Assuming these lumps of gas to be 
compact enough and their number densities small 
enough, so that tidal disruption and collisions are 
avoided over at least a few crossing times, one can 
consider them as super-particles moving through 
the smooth (say NFW) background of DM par- 
ticles. Super-particles experience dynamical fric- 
tion (DF), losing energy to the inner dark halo 
and heating it up. The DF dissipates the energy 
of the gaseous component, but unlike in the case 
of radiative dissipation, the energy is conserved. 
Energy lost by the gas through DF is deposited 
in the inner dark halo and heats it up. For the 
above scenario to work, we have to establish that: 
(a) the coupling via DF is strong enough to pro- 
duce, within a few dynamical times, significant ef- 
fects and (h) with realistic gas mass fraction, there 
is enough energy to heat up the inner halo and 
to form a core, overcoming the competing effect 
of the deepening potential well by the shrinking 
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gas distribution. The semi-analytical Monte Carlo 
model used is described in §2 and the results are 
given in §3. A general discussion of the implica- 
tions and ramifications of the proposed model is 
found in §4. 

2. Method and Model Parameters 

A rigorous modeling of the dynamics of col- 
lapse and virialization of the DM-gas system 
demands hyd-rodynamical/N-body simulations 
starting from cosmological initial conditions and 
with a very large dynamical range, so that the 
DF is properly modeled and the gas dynamics al- 
lows for the fragmentation of the gas component 
into individual massive clumps. Present day ca- 
pabilities, however, only allow for a more modest 
approach. One could, for example, perform high- 
resolution simulations (e.g., Moore et al. 1999; 
Klypin et al. 2000), and introduce ab initio a pop- 
ulation of massive super-particles moving through 
the background DM. Instead, we model the pro- 
posed processes in a semi-analytical Monte-Carlo 
approach, in what amounts to a feasibility study 
of the proposed mechanism. 

Consider a mass M moving with a velocity Vm 
through a homogeneous system of non-interacting 
particles of mass m <^ M and density p with 
a Maxwellian velocity distribution of dispersion 
a. The massive particle experience the DF de- 
celeration (e.g., Chandrasekhar 1943; Binney and 
Tremaine 1987) 




(1) 

where X = Vm / V^cr, the Coulomb logarithm In A 
is defined by A = 6y|f/G(M -|- m) and b is the 
maximal impact parameter of the two body en- 
counters, usually taken to be the radius of the 
system. Although highly simplified and mathe- 
matically incomplete, this purely local treatment 
of DF seems, in many situations of interest, to 
give results which are a reasonable approxima- 
tion to the true behavior (e.g. Zaritsky & White 
1988). The rate of energy loss by a super-particle 
is = M {dV M / dt)-Y M ■ This defines a timescale, 
typically significantly larger than the dynamical 



timescale and related to it through the ratio of 
total mass to t] = Mtot/M. More precisely, at 
any given time X sa 1, allowing one to approxi- 
mate it by rDp/Tdyn ~ (0.75/ln A)?;, where Tdyn 
is the dynamical timescale. In virial equilibrium 
A Ki T], hence TDp/Tdyn ~ 0.75?7/lnr/. The local 
DF timescale will differ significantly from this and 
will be much smaller in the inner dense regions. 

The mathematical model used here is that of 
a spherical halo of mass Mdm with a population 
of Ug super-particles of mass M. While CDM ha- 
los found in N-body simulations are not perfectly 
spherical, they are definitely not flat objects and 
can, for our purposes, be reasonably modeled us- 
ing spherical symmetry. Moreover, the gaseous 
component will be assumed to be in the form of 
compact clumps with sufficiently small crossec- 
tions and extended spatial distribution so that, 
for the timescales considered, dissipative collisions 
leading to the flattening of the gas to a preferred 
plane can be ignored. In spherical symmetry, each 
particle i represents a smoothed shell of the same 
mass rrii and at the same radius r,. Thus, upon 
sorting the particles according to their distance 
from the center, the force on the n^^ particle is 
given by Ti^Z^ Grrii / r'^. To avoid numerical insta- 
bility we use a Plummer softening of 0.1 % the 
initial halo size. 

The numerical model employed here is that of 
Ndm and Ng particles representing the spatial 
distributions of the DM and gas {i.e. super- 
particles), respectively. We choose to determine 
NoM and Ng solely by considerations of a rea- 
sonable statistical representation of a two-fluid 
DM and gas system, coupled via DF. Thus all 
particles have the same mass m = 1/N, where 
N = Ndm + Ng. In the hmit of M > m, the 
dynamical interaction is independent of the value 
of Ndm- To complete the statistical treatment, a 
multiplicative factor F, the 'coupling constant,' is 
introduced in Eq. (1). It determines M and ng in a 
system of total mass Mtot and gas mass fraction fg 
with the same DF interaction strength as our cor- 
responding statistical representation in terms of m 
and Ng. The correspondence is made through the 
following relations: = Ng/N, M = FmMtot 
and Ug = Ng/F. For all the models described 
here the particle numbers are fixed at N = 10^ 
and Ng = 10''. Therefore, fg = 0.1. The code 
is written in units of GMfot = Ij where G is the 
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gravitational constant, so that m = 10~^. The ini- 
tial radius of the halo is taken to be 100 in these 
units. 

The initial DM distribution is constructed to 
follow the NFW density profile p = Psrl/r{rs+r)'^ 
with Vg . Once the total mass is fixed, determining 
rg or ps completely determines the halo model. 
We choose to fix the former quantity (values are 
given in Table 1). The scaling to physical units 
exploits the "universal" scaling relations of NFW 
and is given by 

[r] = 1.63M^/^/i-2/3l0-4 kpc (2) 

and 

[t] = 0.97h-WMDM/Mtot Myr, (3) 

where Mdm is in solar mass and Hubble constant 
is defined by Hq = lOO/i km s~^Mpc~^. Since 
NFW halo density profiles have the same func- 
tional form irrespective of mass, the mass unit 
can be arbitrarily chosen. The system of units 
are completely fixed once the gas mass fraction is. 
One hundred time units will correspond to about 
11 central (in the region r < 10) dynamical times 
and the length unit will be about 1 kpc for a halo 
of 5 X 1O"M0. 

Each Cartesian velocity component is sampled 
from a Gaussian distribution with zero mean and 
variance a = (T{r), where cr is obtained by solv- 
ing the steady state Jeans equation for the ra- 
dial velocity dispersion, the solution of which can 
be obtained in terms of elementary functions and 
quadrature (with boundary condition a = as 
r — > oo). Gas particles are either sampled from 
the same distribution or arc taken to be homoge- 
neously distributed within a certain radius Rg . In 
the latter case, their initial velocities are assumed 
to have a Maxwcllian with half of the halo disper- 
sion at Rg, and the combined system is then scaled 
to virial equilibrium. Systems are left to relax for 
600 to 1000 time units (corresponding to about 6 
to 10 central dynamical times) without including 
the effect of the DF. 

The density and velocity dispersion of the DM 
particles are calculated in a 1000 bins of equal 
numbers and these values are used in evaluat- 
ing the DF acting on the Ng gas particles using 
Eq. (1). The average energy lost by the gas par- 
ticles in a given bin per unit time is updated at 
fixed time steps At = 5. Multiplied by At, this is 



the energy to be gained by the DM particles of the 
same bins. At the end of each time interval, the 
Cartesian velocity components of the DM particles 
arc updated through an additive term chosen from 
a normal distribution with zero mean and variance 
^/ {2/3)Eb/m. Here Eh is the energy gained per 
DM particle. 

The integration is advanced using a Runge- 
Kutta-Merson method with adaptive timestep and 
predetermined local tolcranc;o. The results were 
found to converge for tolerances < 10^'^. Be- 
cause of the strict spherical symmetry, the gravita- 
tional field is much smoother than in the 3D case 
and evolutionary effects are also very slow. After 
reaching dynamical equilibrium (and in the ab- 
sence of frictional forces) , a system of a few thou- 
sand particles may remain virtually unchanged for 
hundreds of dynamical times. The total energy is 
conserved to five digits in the absence of DF. When 
the DF is turned on, the statistical nature of the 
energy feedback and finite size of the shells causes 
additional fluctuations. Still the change in the en- 
ergy did not exceed 0.51% of the total energy, or 
a few percent of the energy exchanged through the 
DF, and was not systematic. 

The results are shown for times up to T = 1800 
(where T = is taken to be the time when the DF 
starts acting) corresponding to ~ 2 Gyr (or about 
18 central dynamical times), although the models 
have been advanced much further. We anticipate 
that eventually the clumps' evolution will become 
dominated by collisions, leading most probably to 
their destruction through star formation, merg- 
ing, fragmentation, and other processes which will 
provide a natural termination to the DF action. 
For example, for Hg objects of radius d randomly 
moving inside a radius Rq, the mean free time is 
l/3rig(-^)^ crossing times. For the parameters of 
Model 1 (c.f. Table ) and assuming d ~ r^/SO for 
example, this is about 25 crossing times. Detailed 
analysis of the aforementioned processes is outside 
the scope of this paper. Here, we either terminate 
the evolution at T = 1800 abruptly (Models 1, 3, 
4), or introduce an exponential cutoff (with char- 
acteristic timescale of 400 units) in the coupling 
parameter F (Model 2). 
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Table 1 
Model Parameters 



Model F Ts Gasatt = Ug M/Mq {ior Mdm = 



1 20 5 uniform Rg = 20 500 2 x 10^ 

2 lOOxe-^/400 5 uniform i?g = 20 100 x e^/^oo 10^ x e-^/^^o 

3 20 5 NFW 500 2 x 10^ 

4 100 3 NFW 100 10^ 



3. Results 

A wide range of models have been calculated 
but only four models are presented here, the pa- 
rameters of which are given in Tabic 1. These 
parameters are representative yet not exhaustive 
of the types of systems we have studied. The gen- 
eral trend is the same for all the models: the DF is 
found to be an efficient mechanism of heating up 
the DM halo, causing the flattening of the density 
cusp into a core-halo structure and consequently 
resulting in less steeply rising rotation curves. In 
principle, it is possible to have a much smaller 
Ng 1000 and to obtain similar results for smaller 

Rgr.5. 

Fig. 1 shows the evolution of quantities charac- 
teristic of the matter distribution of Model 1. We 
note that the final DM rotation curve rises less 
steeply as the DF puffs up the inner halo. To get 
a quantitative assessment of the disappearance of 
the cusp and the emergence of a core-halo struc- 
ture, the final DM rotation curve has been fitted 
by Burkert's (1995) profile, 



^ (r + ro)(r2+r^)' 

claimed to be a good fit to the DM rotation curves 
inferred from observations (e.g., Salucci & Burkert 
2000). 

Assuming the density to be constant, p ~ po 
inside tq and using Eq. 5 of Burkert & Salucci, 
one finds the observational correlation po oc Tq 
between the core radius and density. We now note 
that the change in the rotation curve from a NFW 
profile to a Burkert profile is the result of a de- 
crease in central density accompanied by heating 
of the central region and that the region where 



this mechanism is effective lies inside (as can 
be seen from the lower panels of Fig. 1). Hence, 
one expects ro ~ Vg in the final configuration. In 
addition, for the density to be fit at large radii, 
Po ~ Ps- Moreover, it can be shown (e.g. by 
using the first four equations of Burkert & Silk 
1999 and assuming Ts to be small compared to the 
virial radius) that ps cx r., where 7 = 7.14 

for SCDM models, and 7 = 10 for CDMA cos- 
mologies. Therefore, especially for 7 « 10, one is 
able to anticipate the observational correlation be- 
tween ro and po from the theoretical one between 
Vs and Ps obtained from numerical simulations — 
provided that the DF mechanism is able to wash 
out the density cusps of halos of all masses. 

For comparison, we have studied systems with 
no feedback into the halo from the gas lumps. In 
this expected the halo contracts adiabat- 

ically, becomes more centrally concentrated and 
the kinetic energy profile increases self-similarly 
(Fig. 2). This behavior is thus similar to the case 
of a smooth gaseous system contracting via dissi- 
pational collapse inside the dark matter halo. It is 
then clear that the major assumption here, lead- 
ing to the expansion of the host halo, is that the 
gas is distributed in clumps. The resulting energy 
feedback is able to overcome the competing effect 
of adiabatic contraction. The halo thus expands 
instead of contracting further under the gravita- 
tional influence of the shrinking baryonic compo- 
nent. 

It is possible to obtain nearly identical results 

for larger coupling parameters (thus higher clump 
masses M) acting on much shorter time scales, as 
in the case of Model 2, where an exponential cutoff 
in the two-fluid coupling (parameter) was imple- 
mented. This can be interpreted as representing 
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Fig. 1. — Evolution in central region of Model 1. Top left: Initial and final DM rotation curves. The dashed 
line is a least squares fit using Eq. (4). Top right: Initial and final gas and DM cumulative masses. Bottom 
left: DM density at four different times. Bottom right: Initial and final kinetic energy per DM particle. 



a fragmentation process. In both Models 1 and 2, 
the inflow of clumps due to energy loss via DF is 
small enough so that the total density in the in- 
ner regions also decreases and the total rotation 
curves can still be fit with Burkert models, as can 
be seen from Fig. 3. In Model 1, if our process 
is allowed to continue, it will ultimately cause ex- 
cessive central gas concentration, so that the total 
density fit to Eq. (4) will worsen. This is not true 
of Model 2 however, where due to the exponential 



cutoff, the the coupling via the DF is already too 
small at T = 1800 to have any significant effect. 
In all cases however, the DM distribution contin- 
ues to be fit well by a core-halo profile far beyond 
the timescales shown here. 

The value of F, and hence M, controls the 
rate at which energy is transported, but the in- 
tegral quantity depends mostly on the mass frac- 
tion in clumps in the region of interest (i.e., inside 
r ~ 10). In the above models, after the initial pe- 
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Fig. 2. — Evolution of DM density and kinetic energy per DM particle in Model 1 in the absence of energy 
feedback. 



riod of relaxation, the mass in the gas component 
comprises about a fifth of the mass in this region 
distributed in a much less concentrated way than 
the DM (Fig 1, upper right panel). If the clumps 
are initially distributed the same way as the DM, 
their number density will be too small to have a 
significant effect outside the very central region for 
F = 20. They also have smaller binding energy to 
give. In such cases it was found that the final ro- 



tation curves can be fit by 



^ (r + r,)(^ + r)2' ' 

which rapidly converges to the NFW profile for 
r > Tc — and the requirement that it does so 
fixes the parameter C. Moreover, it was invari- 
ably found that best fits require that rc = A, 
to very high accuracy, making it a one parame- 
ter fit. To compensate for the smaller gas mass 
fraction, one can increase F or the central density 
— which increases the coupling and decreases the 
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Fig. 3. — Top: Initial and final total rotation curves of Model 1 (left panel) and Model 2. Dashed lines 
represent fits using Eq. (4). Bottom: The total density profile for the same models given at four different 
times. 



dynamical time, thus rendering the process much 
faster. In such cases, the final rotation curves at 
t = 1800 can no longer be well fit by Eq. (5): 
the DM distribution is now modified at interme- 
diate radii and is better fit by Eq. (4) (Fig. 4, 
upper panels). Another consequence of the strong 
coupling and condensed initial distribution is that 
the total density profile is very centrally concen- 
trated (Fig 4, lower panels). This is due to the 
more concentrated initial distribution and smaller 



mass, implying smaller amount of binding energy 
that can be released by the clump system, within 
a given radius, before it becomes highly concen- 
trated — it could possibly lead to the formation 
of the galactic bulge/black hole system. This and 
other possible scenarios for the fate of the clumpy 
gas will be discussed elsewhere. 
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Fig. 4. — Top: Initial and final DM rotation curves for Model 3 (left) and Model 4. Dashed lines represent 
best fits using Eq. (5), with i?o = 2.12 and A = 2.14 for Model 3 and i?o = 2.027 and A = 2.017 for Model 4 
(in the units defined in section 2). The dashed dotted line in the right panel represents a fit using Eq. (4). 
Bottom: The total density profiles for the same models given at four different times. 



4. Discussion 

Starting from the NFW cuspy initial density 
profile for the halo, wc have shown that dynamical 
friction provides an efficient mechanism for trans- 
porting energy from the clumpy component to the 
dark matter and heats it up. Under a wide range 
of conditions, the clump orbital energy is sufficient 
to flatten the density cusp and to form a core-halo 
structure. Although we assume that 10% of the 



total mass is in the clumpy component, in prac- 
tice only a fraction of this mass, initially located 
within r ~ rs, participates in the process. This 
is sufficient to change the structure of the central, 
high-density region. 

The general scenario presented here is that of 
the formation of galaxies by the collapse and frag- 
mentation of gas in DM halos. Within the CDM 
paradigm, the DM collapse is followed by the 
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Jeans scale collapse and radiative cooling of the 
gas. The details of this complicated process are 
still largely unknown, but the clumpy nature of 
the gas is unavoidable in the standard cosmogony 
(Hutchings et al. 2000) . It is also required within 
the hierarchical merging scenario of galaxy forma- 
tion. 

For dwarf galaxies, clump masses of 10^ - 10^ 
Mq are required for the DF to be efficient in 
destroying the central cusps. For large galaxies, 
masses of 10^ - 10^ Mq are needed. Such massive 
inhomogeneities are expected to appear after the 
formation of the first objects, when UV radiation 
and inefficient cooling would prevent the subse- 
quent collapse of objects with masses smaller than 
10*^ Mq - 10^ Mq, as argued by Haiman, Rees & 
Loeb (1997). The basic feature assumed here is 
that the gas component is clumpy and that clumps 
are bound enough to survive in the DM back- 
ground for a few dynamical times without colliding 
or disintegrating. If indeed, at least some of these 
massive clumps survive long enough and are not 
tidally stripped or otherwise destroyed, they will 
dissolve the most massive DM cusps. In this way, 
more massive halos would end up with the biggest 
cores — since they initially start with the largest 
NFW scalelength. This is in line with what is in- 
ferred from observations (as shown in section 3). 
Alternatively, a mass spectrum of these clumps 
with a lower mass tail can, in principle, be created 
by a variety of physical processes. The details of 
such processes are outside the scope of this paper 
— although in some of the runs we do allow for 
an exponential cutoff in the gas-DM coupling via 
DF, thus mimicking a fragmentation process. 

If the mechanism described in this paper starts 
acting early in the history of the hierarchical merg- 
ing process, when the halo masses are small, pro- 
portionately smaller clump masses would be suf- 
ficient in changing the CDM halos' central struc- 
tures — since the efficiency of our mechanism de- 
pends on the relative masses of the clumps and 
halo component. In this case the final core radii 
and their variation with halo mass will depend 
on the details of the merging process. The re- 
sultant central structure of these halos, however, 
should still exhibit constant density cores, since 
these are conserved by the merging process (see, 
e.g., Pearce, Thomas & Couchman 1993). 

The model presented here provides a frame- 



work for understanding the possible discrepancy 
between observations of galaxies and the outcome 
of N-body simulations. The final dark halo den- 
sity distributions arc well fit by empirical profiles 
inferred from observations. The less concentrated 
density profiles may also be relevant to the sub- 
structure problem smaller halos with low den- 
sity cores would be less likely to survive to be ob- 
served today. For simplicity, and in order to em- 
phasize the feasibility of the processes described 
here — that the energy deposited in the halo is 
sufficient, given realistic gas mass fraction, in over- 
coming the effect of the adiabatic contraction and 
leading to the formation of a significant constant 
density core — we have considered models with 
zero net angular momentum. Thus we cannot 
comment in detail as to the effect of our mech- 
anism on the angular momentum problem present 
in CDM models of disk formation. In principle 
dynamical friction may also transfer angular mo- 
mentum from the baryonic component to the DM 
one. Our mechanism however is effective only in 
the central region of galaxies and the clumps are 
likely to end up in a relatively low angular mo- 
mentum system, such galactic bulge-central 
mass component. In the outer regions, dynami- 
cal friction is not effective. The fate of the gas 
lumps will be determined by collisions. These will 
lead to collapse to a preferred plane, accompanied 
by contraction and spin up. Once this "second 
generation" material has reached the central re- 
gions it is likely to loose less angular momentum 
to the resulting inner halo — since the latter is less 
concentrated and has higher angular momentum 
— than it would by interaction with the original 
NFW halo. This could lead to less concentrated 
disks. It is also possible that tidal interactions be- 
tween the less concentrated halos lead to higher 
net angular momenta. 

We find that, under certain circumstances, the 
total rotation curves can also exhibit a halo-core 
structure, albeit more modest than the DM alone. 
As emphasized above, one can also envisage a 
scenario where the clumps sink to the bottom 
of the potential well, dissipating their orbital en- 
ergy by the DF and collisions to form the galac- 
tic disk/bulge and the central black hole. Part or 
most of the gas can subsequently be blown away by 
positive energy feedback provided by massive star 
formation and supernovae explosions in the form 
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of galactic winds. The efficiency of this process 
should anti-correlate with the depth of the poten- 
tial well and hence distinguish between dwarf and 
massive galaxies. 

We use a semi-analytical Monte-Carlo ap- 
proach, based on the Chandrasekhar approxi- 
mation for the DF drag force in a homogeneous 
medium of uncorrelated particles. This approxi- 
mation can be improved by high-resolution simu- 
lations in which the super-particles are introduced 
by construction and where the full iV-body inter- 
action, including potentially important global re- 
sponse of the background medium (e.g., Weinberg 
2000) is taken into account. N-body simulations 
of the dynamics of satellite galaxies in massive 
halos show that the semi-analytical modeling un- 
derestimates the DF effect when the satellite mass 
fraction is 0.05 — 0.2 (Tormen et al. 1998), similar 
to the fraction of baryonic mass within r < 10 
in our models. The eventual slowing down of the 
self-consistent evolution compared to the semi- 
analytic one in their models is probably due to 
the decreasing central background density result- 
ing from energy feedback into the background 
medium from the dissipating satellites, which is 
not modeled in their semi-analytic calculations. 
We have also assumed spherical symmetry. Sim- 
ulated DM halos arc found to be triaxial. In 
such systems the effect of DF is strongly enhanced 
(Pesce, Capuzzo-Dolcetta & Vietri 1992). There- 
fore, the effect described in this paper can be 
probably achieved for substantially smaller clump 
masses or on smaller timescales. 
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